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Context. It is commonly believed that the observed overabundances of many chemical species relative to the expected 
cluster metallicity in blue horizontal branch (BHB) stars appear as a result of atomic diffusion in the photosphere. The 
slow rotation of BHB stars (with T e g > 11,500K), typically wsini < 10 km s _1 , is consistent with this idea. 
Aims. In this work we search for observational evidence of vertical chemical stratification in the atmosphere of HD 135485. 
If this evidence exists, it will demonstrate the importance of atomic diffusion processes in the atmospheres of BHB 
stars. 

Methods. We undertake an extensive abundance stratification analysis of the atmosphere of HD 135485, based on 
recently acquired high resolution and S/N CFHT ESPaDOnS spectra and a McDonald-CE spectrum. 
Results. Our numerical simulations show that nitrogen and sulfur reveal signatures of vertical abundance stratification 
in the stellar atmosphere. It appears that the abundances of these elements increase toward the upper atmosphere. This 
fact cannot be explained by the influence of microturbulent velocity, because oxygen, carbon, neon, argon, titanium 
and chromium do not show similar behavior and their abundances remain constant throughout the atmosphere. It 
seems that the iron abundance may increase marginally toward the lower atmosphere. This is the first demonstration 
of vertical abundance stratification of metals in a BHB star. 
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1. Introduction flow onto HD 135485 from a high mass component may 

have created the observed enhancements of metal species. 
HD 135485 appears to be an evolved sta r that lies close to or In order to try to prove this hypothesis, Dufton investigated 
5, ■ on the horizontal branch (Trundle et al.EOOl). It has mid-B the radial velocity during a nve month period, but found 
03 ; spectral type and shows many prominent absorption lines n0 variation (see Trundle et al. (j2TKJll for details). 

of metallic species in its spectrum. The most recent abso- « , , . . , r 

. , . ..a? ,. , . f . . r TT „ , or . nr Another possible explanation proposed tor the unusu- 

lute and differential abundance analysis ol HD 135485 was , • , , ,,. ., . ,f , , c , • 

nerformed bv Trundle et al (rMUl Thev found a laree (0 5 ally hlgh metallicit y 1S that the star was formed m a re S 10n 
pciioiiiieu uy n miuic ci ai. izuuii. r ncy iouuq a iait;c iu.o- ri • 1 , tv 1 j_i /~i 1 . * r-i .l j j 

* . . , / , . *—~ i — ' , . ° \ ol high metallicity, such as the Galactic Center, and during 

1.0 dex enhancement ol metal abundances m comparison , . , 7... , , ,. , , j 

, ... „. ., „ 1, • 1 the long period ol its stellar evolution it could have traveled 

with the solar composition. Similar results were obtained , j- , / ocl \ 1 , , 1 1 rp 11 

. . „ . .. , . _ . r over a large distance (^26 kpc as demonstrated by Irundle 

previously by Schonberner and Dufton from ^ &} mm from ^ . g of ^ kinematical propert ies of 

analysis of moderate-resolution photographic spectra. In an ^ 35485 

attempt to explain the anomalous chemical composition in 

the atmosphere of HD 135485, Dufton (1575)1 considered Together with general enrichment of all metals (except 
the possibility of binary system evolution, where the mass mckel which 18 significantly depleted) with respect to hy- 
drogen, Trundle et al. ()2001|) have also found a slightly en- 
hanced helium abundance of ~0.3 dex. They interpret this 



Send offprint requests to: V. Khalack fact along with the anomalously high (~1.3 dex) nitrogen 

e-mail: khalakvOumoncton. ca overabundance as an indication that the carbon-nitrogen 

Based on observations made with ESPaDOnS at the (CN) le took lace duri the hydrogen-burning phase 

Canada- France-Hawaii Telescope (CFHT) operated by the A m, , . . , , , 
„ ,. , „ , „ ., ,, m Vn v „ „ 1 ii T j. 1 on the mam sequence. Ihe relatively normal carbon cn- 
JNational Research Council (NRC) of Canada, the Institut des . -in ^ 1 1 , 1 ■ 1 
Sciences de l'Univers of the Centre National de la Recherche hancement might reflect a balance between its being de- 
Scientifique (CNRS) and the University of Hawaii and on ob- stroyed by the CN cycle and created via helium core 
servations made with Echelle Spectrograph on the McDonald flashes during the red giant phase (which could also ex- 
Observatory 2.1-meter Otto Struve Telescope. plain the observed neon enrichment in the atmosphere). 
** Current address: US Naval Observatory, 3450 Massachusetts These facts and the present position of HD 135485 on the 
Avenue NW, Washington DC 20392 Hertzsprung-Russell (HR) diagram argue that this star is 
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Table 1. Journal of spectroscopic observations of HD 135485. 



Date 


HJD 


Instrument 


^cxp 


S/N ratio 


Coverage 


Resolution 


(UT) 


(2450 000 +) 




(s) 


(pixel)" 1 


A (A) 




2001 Mar 09 


1978.9242 


CE 


600 


67 


4760 - 5570 


60000 


2005 May 19 


3509.8677 


ESPaDOnS 


600 


150 


3720 - 10290 


80000 


2005 May 19 


3509.8753 


ESPaDOnS 


600 


140 


3720 - 10290 


80000 



evolved and probably belongs to the blue horizontal branch 
(BHB) stars. The observed low (projected) rotational veloc- 
ity (v s'mi < 4 km s _1 ) of HD 135485, the observed chemical 
peculiarities (except the helium enrichment) and the effec- 
tive temperature T cff = 15500K (Trundle et al. I200ip are 
typical for BHB stars. 

Comprehensive surveys of BHB star abundances show 
that the stars hotter than T e g ~ 11,500K have abundance 
anomalies when compared to the BHB stars of lower effec- 
tive temperature in the same clusteiQ (Glaspey et al.[l989, 
Grundhal et al. I1999j) . Also, Peterson, Rood & Crocker 
(fT995|) ; Behr et al. (|2000)) and Recio-Blanco et al. (j2004|) 
showed the existence of a discontinuity in stellar rotation 
velocity distribution of BHB stars at T cff ~ 11,500K. All 
of the hotter stars show modest rotation (usinz < 10 km 
s -1 ), while the cooler stars are rotating more rapidly. 

Other observational results show both photometric 
jumps (Grundhal et al. I1999|) and photometric gaps 
(Ferraro et al. I1998p . also occurring at T e g ~ 11,500K, 
as well as low measured surface gravities. The low gravity, 
abundance anomalies and slow rotation suggest that micro- 
scopic atomic diffusion is effective in stellar atmospheres of 
the BHB stars with T c g > 11,500K. In this picture, the com- 
petition between radiative acceleration and gravitational 
settling yields a net acceleration on atoms, which results 
in their diffusion through the atmosphere. A consequence 
of this process is vertically-stratified abundances of chemi- 
cal elements. The aim of this paper is to attempt to detect 
obscrvationally such vertical stratification of elements in 
the atmospheres of BHB stars. Direct measurement of this 
stratification from line profile analysis would provide a con- 
vincing argument in favour of efficient chemical diffusion in 
the atmospheres of BHB stars with T e ff > 11,500K. 

The most important difference between HD 135485 and 
most BHB stars is that it has an enhanced helium abun- 
dance. Therefore HD 135485 might be a very young BHB 
star in which helium has not yet settled. However, a de- 
tailed investigation of this possibility is outside the scope 
of this paper. 

To verify if the abundances of some chemical species are 
vertically stratified in the atmosphere of HD 135485, we 
undertook a detailed analysis of available high-resolution 
spectra employing a modification of the approach for abun- 
dance stratification analysis developed by Ryabchikova et 
al. ([2005). In Sect. [2] we discuss the properties of the ob- 
served spectra, and in Sect. [3] we describe the atmospheric 
parameters of HD 135485 and modelling details. The re- 
sults of abundance determinations are considered in Sect.lH 
while a description of vertical stratification for some chem- 
ical species is shown in Sect. [5] A discussion of our results 
follows in Sect. [6] 



2. Observations 

Spectroscopic observations of HD 135485 were undertaken 
in March 2001 and May 2005. The journal of spectroscopic 
observations is shown in Table [TJ where individual columns 
give the UT date of the observation, heliocentric Julian 
Date of the observation, the instrument specification, expo- 
sure time, signal-to-noise (S/N) ratios per pixel (the peak 
typically occurs around 5400 A), spectral coverage and re- 
solving power (which depends on central wavelength) . 

2.1. McDonald Cassegrain Echelle spectrum 

The first spectrum of HD 135485 (from Table HJ was 
obtained in the framework of a study of the rotational 
characteristics (wsini) of BHB stars (Behr I2003|) using 
the Cassegrain Echelle (CE) Spectrograph (McCarthy et 
al. [T593) on the McDonald Observatory 2.1-meter Otto 
Struve Telescope. This instrument provides a nominal re- 
solving power of R ~ 60000. The signal-to-noise ratio (per 
pixel) was computed from the reduced spectrum by smooth- 
ing the spectral regions identified as continuum, and then 
computing the rms deviation between the observed and 
smoothed spectra. 

The package of routines developed by McCarthy (|1990p 
for the FIGARO data analysis package (Shortridge ]1993) 
was employed to reduce the spectrum. The detailed descrip- 
tion of the reduction procedure for this particular spec- 
trum, including the wavelength calibration routine and the 
continuum normalization, is provided by Behr (2003). The 
procedure follows a standard prescription, with bias sub- 
traction, flat-fielding, order extraction, and wavelength cal- 
ibration from thorium-argon arc lamp observations. Cosmic 
ray hits were identified and removed by hand for maximum 
spectral quality. No smoothing was applied to the spectra. 

2.2. CFHT-ESPaDOnS spectra 

The other two spectra were obtained using the new 
ESPaDOnS (Echelle SpectroPolarimctric Device for 
Observations of Stars) spectropolarimeter at the Canada- 
France-Hawaii Telescope (CFHT). ESPaDOnS is fun- 
damentally similar in construction to the MuSiCoS 
spectropolarimeter (Donati et al. 1999) and allows the 
acquisition of an essentially continuous spectrum through- 
out the spectral range 3700 A to 10500 A in a single 
exposure. The optical characteristics of the spectrograph 
as well as the spectropolarimeter observing procedures, are 
described by Donati et al. (in preparation^- Observations 
were performed in the spectroscopic 'sky only' mode of 
the instrument that provides the highest resolving power 



1 Most of the known BHB stars are found in globular clusters. 



2 For more details about this instrument, the reader is invited 
to visit www.ast.obs-mip.fr/projects/espadons/espadons.html 
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(R ~ 80000). This resolution corresponds to velocity 
elements of 3.75 km s _1 . 

The ESPaDOnS spectra were reduced using the Libre- 
ESpRIT reduction tool (Donati et al., in preparation), 
which is the most recent release of ESpRIT (Donati et al., 
1997). The continuum of each spectral order was normal- 
ized by fitting a 3rd to 5th order polynomial (as selected 
by eye) to regions free of prominent spectral lines and not 
located immediately at the beginning or at the end of the 
spectral order. 

3. Line profile simulations 

3.1. Stellar atmosphere parameters 

The line profile simulations were performed using a Phoenix 
(Hauschildt et al. [TWTj) LTE (Local Thermodynamic 
Equilibrium) stellar atmosphere model with enhanced 
(+1.0 dex) metallicity. The model is calculated assum- 
ing T c g = 15, 500K, and logg = 4.0, which according to 
Trundle et al. (|2001[) provide the best description of ob- 
servable line profiles of HD 135485. Those authors adopted 
a microturbulent velocity of 5 km s _1 to provide a con- 
stant behavior of abundances over the range of equivalent 
widths measured for sulfur and nitrogen ions. Behr (2003) 
has adopted for this star £ = 1.4 km s _1 and obtained 
■usini = km s _1 . 

From our perspective, the physical motivation for 
adopting a microturbulent velocity is unclear. According 
to Landstreet (1998), photospheric convection should be 
very weak at the effective temperatures of BHB stars. 
Landstreet found that spectroscopically-derived microtubu- 
lence had decreased from around 4 km s _1 in late A stars 
with T cff ~ 8000 K to 0.5 km s _1 in early A stars with 
T e ff ~ 10500 K. In particular, strong microtubulence (e.g. 
5 km s" 1 as adopted by Trundle et al. 2001) should be ac- 
companied by strong line asymmetries, which we do not ob- 
serve. We therefore realistically expect no detectable micro- 
turbulence in HD 135485. Of course, Landstreet's analysis 
was performed for main sequence stars; we might reason- 
ably expect that BHB stars show some differences, although 
the atmospheric parameters of HD 135485 are very similar 
to those of main sequence stars. 

Based on this guidance, we have tentatively adopted a 
microturbulent velocity km s _1 for our abundance anal- 
ysis of HD 135485. Using zero microturbulence in spectral 
synthesis we estimated stellar rotational velocity vsini =2 
- 3 km s^ 1 . This is an important assumption for several rea- 
sons. First, because diffusion is a slow process, the mixing 
which might be implied by the presence of microturbulence 
could serve to reduce or erase chemical stratification. Also, 
adoption of microturbulence modifies line strengths, lead- 
ing to systematic changes in inferred mean abundances as 
well as abundance trends with line strength. As we shall 
see in Sect. 5, the validity of the assumption of zero micro- 
turbulence can be tested, and it is fully consistent with our 
spectroscopic observations. 

3.2. Procedure 

All of the available spectra were examined to compose a 
line list for chemical species that are suitable for abun- 
dance and stratification analysis. The line identification 
is performed using the VALD-2 (Kupka et al. 1999: 



Table 2. List of spectral lines used for the abundance anal- 
ysis. 



Nil 


A, A 


log gf 


Ei,cm 1 


log 7rad 


Ref. 


4447.030 


0.228 


164610.76 


9.16 


NIST 


4507.556 


-0.817 


166678.64 


9.33 


VALD 


4601.478 


-0.428 


148940.17 


9.23 


VALD 


4613.868 


-0.665 


148940.17 


9.23 


VALD 


4621.396 


-0.514 


148940.17 


9.23 


VALD 


4630.539 


0.094 


149076.52 


9.15 


VALD 


4779.722 


-0.587 


166521.69 


9.58 


VALD 


4803.287 


-0.113 


166678.64 


9.58 


VALD 


4987.376 


-0.555 


168892.21 


9.34 


VALD 


4994.367 


-0.069 


168892.21 


9.34 


VALD 



Ryabchikova et al. fT999l and NISlEI (version 3.0.3) line 
databases (Reader et al. I2002|) . The same sources are used 
to specify the atomic data for the selected lines. For Fell 
and Cm lines we extracted atomic data from Raassen & 
Uyling^3 (|1998j) , while for some Fe n lines we also used the 
atomic data recently obtained by Fuhr & Wiese l|2006p . An 
example of the first ten lines from the line list selected for 
the abundance analysis of N n lines is shown in Table Iffil 
where references identify sources for the adopted log gf val- 
ues, while the given energy levels are taken mostly from the 
NIST database. 

All line profiles are simulated using the Zeeman2 spec- 
trum synthesis code (Landstreet 119881 Wade et al. [2001 ) 
assuming Gaussian instrumental profiles. The stability of 
the ESPaDOnS instrumental profile has been investigated 
using the calibration spectra of a ThAr quartz lamp. We 
find the instrumental profile to be stable during the night 
of observation and to reveal no dependence on the inten- 
sity of the line profile. The Zeeman2 code has been mod- 
ified (Khalack & Wade I2006|) to allow for an automatic 
minimization of the model parameters using the downhill 
simplex method (Press et al. I1992[) . The relatively poor ef- 
ficiency of the downhill simplex method, requiring a large 
number of function evaluations, is a well known problem. 
Repeating the minimization routine several times in the 
vicinity of a supposed minimum in the parameter space 
allows us to check if the method converges to the global 
minimum. 

To search for the presence of abundance stratification, 
we have used two different methods. In the first method, 
we estimate the abundance of a chemical element from in- 
dependent analysis of each line profile. For every layer of 
the stellar atmosphere model (it contains 50 layers) we cal- 
culate the line optical depth T£ in the line core. We sup- 
pose that each analyzed profile is formed mainly at Tt=l, 
which corresponds to a particular layer of the stellar at- 
mosphere. Finally, for this layer the respective continuum 
optical depth T5000 is specified. In this way, we have built 
the scale of optical depths T5000 to track vertical abundance 
stratification through the analysis of all the available line 
profiles for a chosen chemical element. 



3 http:/ /physics. nist.gov/PhysRefData/ASD/index. html 

4 |ftp:/ / ftp.wins.uva.nl/pub /orth| 

J Table [2] is presented full in the electronic edition of 
Astronomy and Astrophysics. A portion is shown here for guid- 
ance and content. 
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Three free model parameters (the element's abundance, 
line radial velocity V r and wsini) were derived from each 
line profile using the aforementioned automatic minimiza- 
tion routine. HD 135485 has very sharp absorption lines 
and even for resolution R=80000 a line profile usually con- 
sists of only 8 to 15 bins. This leads to the fact that the 
downhill simplex method fails to converge to the global 
minimum in approximately 2 to 5% of analyzed line pro- 
files. This means that the automatic minimization routine 
remains in the vicinity of a local minimum and results in a 
higher final value of the x 2 -function. Repeating the mini- 
mization routine with different starting values of the model 
parameters usually solves this problem and provides a bet- 
ter fit, which is checked using a few additional runs of the 
simulation routine. Sometimes, analysis of a line profile re- 
sults in a radial velocity which differs significantly from the 
average V T =4. 2 km s -1 (obtained from analysis of all line 
profiles of the same chemical element and other chemical 
species). This may be evidence for line misidentifications 
or inaccurate line wavelengths. Therefore lines with radial 
velocities which differ more than 1.5 km s _1 from the aver- 
age are not taken into account in the modeling and are not 
included in Tabled 

In the second method, we only analyzed those lines 
which appear to be good candidates in the previous ap- 
proach and consider the two-step stratified abundance 
model with a linear transition zone (Ryabchikova et al. 
2003 ; Wade et al. 2003) to describe the abundance distribu- 
tion. In this case, all the selected lines were analyzed simul- 
taneously and in the minimization routine we operated with 
6 free model parameters: the element's abundance in the 
deeper layers, the (standard) optical depth n of the tran- 
sition zone lower (deeper) boundary (where the element's 
abundance begins its linear increase or decrease), the op- 
tical depth T2 of the transition zone upper boundary, the 
difference between the abundances at two boundaries, the 
stellar radial velocity and the vsmi. We note that for the 
second method the downhill simplex method never failed 
to reach the global minimum. This means that information 
stored in all of the analyzed lines, modeled simultaneously, 
is statistically sufficient to provide a unique solution. To 
determine the averaged abundances for all of the analyzed 
chemical species, we simultaneously fit the available line 
profiles using the same three free model parameters used in 
the first approach. 

Initially, we specify arbitrarily the free model parame- 
ters in the range of their possible values and simulate the 
line profiles. We then compare the simulated profiles with 
the observed spectra, and calculate the reduced x 2 ■> which 
we adopt as a measure of the fit quality. The expression 
for the \ 2 function, which reflects the agreement between 
the simulated profiles (ii,Aj) and observed spectra (I"^) is 
given by: 




where cr[I° b ^] corresponds to the measurement errors, Nj 
represents the number of profiles indexed by i, while Ni 
is the number of pixels in each analyzed line profile. In 
the first method of the analysis we used Nj=l. In fact, the 
simulated spectra are calculated with a resolution of 0.01 A 
and so do not provide a direct coincidence of wavelengths in 



Table 3. Mean atmospheric abundances log(iV/iVtot) °f 
HD 135485 with uncertainties equal to the standard de- 
viation of n measured lines. 



Ion 


CFHT-ESPaDOnS 

log(JV/JVtot) n 


McDonald-CE 

log(JV/JV t ot) n 


Sun 

log(JV/JV, 


He I 


-0.58±0.16 


9 


-0.54±0.15 


3 


-1.11 


C. I 


-2.78±0.17 


10 






-3.65 


Cn 


-2.69±0.14 


15 


-2.88±0.16 


4 


-3.65 


Ni 


-2.60±0.20 


3 






-4.26 


Nil 


-3.00±0.17 


21 


-3.01±0.18 


9 


-4.26 


Oi 


-2.82±0.13 


11 


-2.74±0.19 


2 


-3.38 


On 


-3.28±0.22 


9 






-3.38 


Nei 


-3.22±0.06 


4 






-4.20 


IN a I 


-4.22±0.24 


2 






-0.8 I 


Mgi 


-3.62±0.20 


1 






-4.01 


Mgn 


-3.73±0.12 


4 






-4.51 


Alii 


-5.32±0.11 


7 






-5.67 


Aim 


-5.10±0.09 


4 






-5.67 


o; tt 
31 II 


-4.28±0.23 


9 


-4.23±0.11 


2 


-4.03 


Sim 


-3.96±0.22 


2 






-4.53 


Pn 


-5.81±0.12 


8 


-5.76±0.13 


3 


-6.68 


Sn 


-4.43±0.18 


48 


-4.47±0.17 


21 


-4.90 


CI ii 


-6.35±0.20 


1 


-6.47±0.20 


1 


-6.54 


Aril 


-5.03±0.05 


4 


-4.96±0.19 


2 


-5.86 


Tin 


-6.70±0.04 


5 






-7.14 


Cm 


-6.01±0.07 


5 


-5.89±0.14 


2 


-6.40 


Fen 


-4.08±0.16 


63 


-4.16±0.18 


31 


-4.59 


Srii 


-8.50±0.20 


1 






-9.12 



the simulated and observed spectra. Therefore, during the 
X 2 function evaluation the simulated spectral intensity at 
the exact observed wavelength is calculated using a linear 
interpolation. 

In order to evaluate the fit errors (and therefore the 
uncertainties on the derived free parameters), we calculate 
deviations of the simulated profiles produced as a result of 
small variations of each of the free parameters, thus intro- 
ducing a small shift along one axis in the x 2 hyper-space 
from the point of the function minimum value. Using this 
procedure, and taking into account the uncertainties of the 
observational data and the obtained minimum value of the 
X 2 -function, we can estimate the errors of the best-fit pa- 
rameters. 

4. Results of abundance analysis 

We have identified in the spectra of HD 135485 most of 
the lines reported by Trundle et al. ((2001). No lines of 
Mn II and Ni II were detectable in our spectra and these ele- 
ments were therefore not analysed. It seems that Trundle et 
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al. (|200ip have determined the Mnil and Nin abundances 
from the much stronger ultraviolet lines of these species. 
The McDonald-CE spectrum has lower spectral resolution 
than the ESPaDOnS data and is used mainly to confirm the 
average abundance obtained with the ESPaDOnS spectra. 

A preliminary LTE abundance analysis of helium lines 
suggests an enhanced helium abundance in HD 135485, and 
the element also seems to be vertically stratified. The he- 
lium overabundance obtained is in good agreement with the 
results reported by Trundle et al. (|200ip . 

The mean photospheric abundances we derive for 
HD 135485 are reported in Tabled The second and third 
columns contain respectively the mean abundance of the 
chemical species and the number of analyzed line pro- 
files in the ESPaDOnS spectra, while the fourth and fifth 
columns represent similar results obtained from analysis of 
the McDonald-CE spectrum. The last column contains the 
solar atmospheric abundances recalculated from \og(N / Nn) 
data (Asplund et al. lMJS]) . 

4.1. Light elements: C to Ar 

Carbon is enhanced by ~0.9 dex in comparison with the so- 
lar abundance and does not show strong evidence of vertical 
abundance stratification (see Fig. UK). Nitrogen appears to 
be strongly enhanced by ~1.7 dex for the neutral species 
and by ^1.3 dex for the first ion. Unlike C, it shows sig- 
natures of vertical abundance stratification (see Fig. [2^). 
The abundances inferred from singly-ionized oxygen is al- 
most the same as its solar value, while the neutral species 
is enhanced by ~0.5 dex. 

Neon and argon appear to be enhanced respectively by 
~1.0 dex and ~0.8 dex in comparison with their solar abun- 
dances. The sodium abundance is estimated from the anal- 
ysis of two comparatively weak Nai lines at 5688.205 A and 
8194.824 A, which results in a ~1.6 dex enhancement. 

Only one unblended Mgi line (5172.684 A) was found, 
which results in a magnesium enhancement of ^0.9 dex. 
Meanwhile, singly ionized magnesium is enhanced by ^0.8 
dex. Alii and Aim are slightly enhanced by ~0.3 dex and 
~0.5 dex respectively. Singly ionized silicon has only a ~0.2 
dex enhancement in comparison with its solar abundance, 
while Si in is enhanced by ~0.6 dex. Phosphorus shows ~0.9 
dex enhancement. Sulfur appears to be vertically strati- 
fied like nitrogen (see Fig. [^b) and its average atmospheric 
abundance is enhanced by ~0.4 dex. Chlorine seems to have 
a solar abundance. 

4.2. Iron-group elements: Ti to Fe 

The spectrum of HD 135485 is extremely rich in Fell 
lines. Many of them originate from high-excitation levels 
with Ej > 10 eV. For Fen lines, we have used mainly 
Raassen & Uylings (|1998j) results to specify their oscil- 
lator strengths and the VALD-2 database to specify the 
damping coefficients. The Raassen & Uylings (1998]) os- 
cillator strength data usually provides a slightly better 
agreement between observed and simulated Fen line pro- 
files than the Kurucz (g9 93|) GF IRON list (Pickering et al. 
I200H Ryabchikova et al. [2005). The obtained average at- 
mospheric Fell abundance is enhanced by ~0.5 dex and 
shows some signs of stratification (see. Fig. [T|d), but in 



an opposite manner to nitrogen and sulfur (i.e. with the 
concentration increasing toward deeper layers). 

Using the ESPaDOnS spectra we find that chromium 
and titanium are represented by only a few unblended lines 
and both elements appear to be enhanced by ~0.4 dex. 
The oscillator strengths for Cr II lines have been taken from 
Raassen & Uylings ()1998[) . while for Ti n lines we have used 
VALD-2 database. 

We also checked the Sr II line at 4077.71 A, which is only 
marginally visible in these spectra. Its abundance analy- 
sis reveals a ^0.6 dex enhancement - much less than the 
value 1.24 dex enhancement obtained by Trundle et al. 
(2001). Nevertheless, the enhancements of the other chem- 
ical species reported in Table [3] are in good agreement with 
the results of Trundle et al. (|200Tj) . 

Adoption of different values of microturbulent velocity 
by Trundle et al. (2001 ) and by us should lead to systematic 
differences in our derived abundance. However, these differ- 
ences may be reduced because we have used different line 
lists and atomic data for the abundance analysis. Trundle 
et al. (|200ip have used mostly NIST atomic data for all 
chemical species, while we have tried to use the most pre- 
cise atomic data currently available for each element (see 
Table [2]). 

5. Vertical abundance stratification 

Applying the technique described in Sect. 13.21 we have tried 
to determine if the abundances of some chemical species are 
vertically stratified. To reach this aim, we need to analyze as 
many line profiles with different depths of formation (char- 
acterised by different excitation potential and line strength) 
as possible. According to Table [3J a comprehensive anal- 
ysis can only be performed for carbon, nitrogen, oxygen, 
sulfur and iron, which are represented in the spectrum of 
HD 135485 by a sufficient number of lines. In general, we 
have selected for our analysis lines free of predicted or in- 
ferred blends. However, if a blend is from a line of the same 
chemical element that forms the main line profile, such a 
line was also included in our simulation. 

The very sharp absorption lines of the HD 135485 (8 to 
15 bins per line profile), the relatively low signal-to- noise 
ratio and the uncertainties in the atomic data provide com- 
paratively high errors in abundances inferred from a single 
line. Taking into account these uncertainties as well as ion- 
ization balance errors, we can still see some tendencies of 
vertical abundance stratification for some of the analyzed 
chemical elements. Fig.[5]shows a systematic trend, wherein 
N and S lines formed higher in the atmosphere provide 
abundances which are significantly larger than those formed 
lower in the atmosphere. To check the influence of differ- 
ences in log g/ values in different databases we have per- 
formed independent simulations of sulfur lines using NIST 
and VALD2 atomic data and have found similar systematic 
trends in both cases. 

The trend observed in the abundances of N and S was 
previously reported by Trundle et al. (2001), but they in- 
terpreted it in terms of microturbulence. However, such an 
interpretation would require a similar behaviour of the lines 
of C, O and Fe, which we do not observe. For example, al- 
though the inferred nitrogen and sulfur abundances appear 
to increase systematically with line strength, carbon and 
oxygen abundances are approximately constant, whereas 
the iron abundance shows a weak tendency to decrease with 
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Fig. 1. Abundance estimates from the analysis of a) Ci (filled circles) and Cn (open circles) and b) Fen (open circles) 
line profiles as a function of line (core) formation optical depth. No clear signatures of abundance stratification are visible. 
The abundances obtained have approximately the same value for all optical depths at which carbon lines are formed. 

a) b) 
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Fig. 2. Same as in Fig. [T] but for a) nitrogen and b) singly ionized sulfur. Ni (filled circles) and Nil (open circles) lines 
which are formed in layers with log T5000 <-1.0 require an enhanced abundance to be fit in comparison with lines of the 
same ionization stage formed deeper in the atmosphere. A similar situation is found for S 11 (open circles) lines formed 
in layers with log T5000 <-1.5. 



Table 4. Parameters of the vertical abundance distribution. Here the columns specify respectively the element, its 
abundance in deeper layers, the abundance difference between the two boundaries, the optical depth t\ of the transition 
zone's lower (deeper) boundary, where the element's abundance increases (or decreases) linearly, the optical depth T2 of 
the transition zone's upper boundary, the radial velocity and v sinz with the respective error estimates and the fit quality 
for the model with abundance stratification Xs an d without it Xn- 



El. 


Abundance 


in log(JV/JVtot) 


Atmospheric 


depths in log t 50 oo 


Vr 


V sini 








Low atmosphere 


Abundance difference 


logTl 


logT2 


km s _1 


km s _1 


2 

Xs 


2 

Xn 


N 


-3.15±0.03 


0.68±0.25 


-0.67±0.27 


-2.75±0.45 


-4.20±0.13 


1.3±0.7 


2.61 


3.83 


S 


-4.63±0.02 


0.94±0.28 


-1.31±0.16 


-3.00±0.40 


-4.24±0.08 


3.0±1.2 


6.87 


10.06 



line strength (see Fig.[T]). The lack of a coherent trend in the 
abundances derived from weak versus strong lines in spectra 
of different elements is not consistent with a microturbulent 
interpretation, and is fully consistent with chemical strati- 
fication. Our simulations of Fe 11 lines show that adoption 
of a microturbulent velocity £ = 2 km s _1 results in well 



visible decreasing of iron abundance with line strength. We 
point out that this conclusion does not necessarily imply 
that microturbulence is not present, but rather that it can- 
not explain the observed systematic abundance trends, that 
its presence is not implied by any of the data discussed in 
this paper, and that its presence would be surprising given 
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published observational and theoretical results (Landstreet 
1998). Our results therefore suggest that the observed sys- 
tematic trends in the abundances of N and S imply that 
these elements are vertically stratified in the atmosphere of 
HD 135485. 

To obtain an estimate of the abundance distribution 
as a function of optical depth for N and S, we have em- 
ployed a two-zone empirical model with a transition zone 
(Ryabchikova et al. I2U031 Wade et al. 12003]) and analyzed 
simultaneously the same list of lines selected for each ele- 
ment. The final results for nitrogen and sulfur stratification 
in the atmosphere of HD 135485 are presented in Table [4] 
and are illustrated in Fig. [3] They are in a good accor- 
dance with the abundance stratification data that we have 
obtained from the first method for N and S (see Fig. [2]). For 
comparison, we have also simulated the same list of lines 
assuming the same abundance of the element (fitted as a 
free parameter) at different atmospheric depths. From the 
last two columns in Tabled] we can see that the model with 
stratified abundance distribution results in a 30% lower \ 2 - 
function than the model with uniform abundance. This fact 
shows that application of stratified model significantly im- 
proves our fit of analysed line profiles. 

The N and S abundance stratifications were also in- 
ferred assuming other stellar atmosphere models with the 
same gravity log g = 4.0, but with the different effective 
temperatures T e g = 14, 500K and T c ff = 16, 500K, to es- 
timate the influence of effective temperature errors on the 
vertical abundance stratification profiles. Our simulations 
show that an effective temperature increase does not ef- 
fect the final results, while its decrease leads to significant 
changes of the vertical abundance stratification profile. The 
stellar atmosphere model with lower temperature results 
in a steeper stratification profile for S and an increase in 
the average S abundance by 0.4 dex, while for nitrogen it 
provides approximately a 0.8 dex increase in abundances 
derived from lines formed deeper in the atmosphere where 
logTsooo > —1-0. This leads to a reversed stratification pro- 
file for N. Therefore the results obtained (see Table [4]) for 
sulfur can be considered to be a lower limit for its abun- 
dance stratification profile, whereas a uniform vertical dis- 
tribution of N cannot be ruled out. 

In conclusion, the analysis of our spectra shows vertical 
stratification of of S and possibly N, while the elements 
C, O, Ne, Ar, Cr and Ti seem to have an homogeneous 
abundance throughout the parts of the atmosphere that we 
can sample. Meanwhile, we cannot draw firm conclusions 
concerning the stratification of Fe, Mg, Si, P and Al. Some 
of these elements are represented by a small number (see 
Table [3]) of unblended line profiles and it is hard to obtain 
clear results. 

6. Discussion 

According to Trundle et al. (]2"DUT)) . HD 135485 is an evolved 
star that probably belongs to the BHB group. Its slow ro- 
tation, high effective temperature (T e g=15,500 K) and ob- 
served abundance anomalies (see Table [3]) support the idea 
that microscopic atomic diffusion is effective in its atmo- 
sphere. 

Our analysis suggests that sulfur is moderately over- 
abundant and increases in concentration in the upper at- 
mospheric layers. Nitrogen is strongly overabundant and 
shows a similar signature of vertical stratification, although 




-3 -2.5 -2 -1.5 -1 -0.5 

logTsooo 



Fig. 3. Abundance stratification in the atmosphere of 
HD 135485 for nitrogen and sulfur. 

the derived stratification of this element depends sensitively 
on the adopted effective temperature (which is uncertain by 
about ±1000 K). Similar results for these two species were 
obtained earlier by Trundle et al. (|200ip . but explained in 
terms of a microturbulent velocity field. Absence of any 
abundance stratification for several other elements as well 
as some hints that the iron abundance might be enhanced 
in deeper atmospheric layers (i.e. opposite to the tendency 
shown by N and S) shows that this interpretation is not 
able to correctly explain the spectroscopic properties of this 
star. A more plausible explanation is that the abundances 
of N and S vary as a function of depth within the atmo- 
sphere. The detection of a similar trend of Fe abundance 
versus depth in the BHB star WF4-3085 in the globular 
cluster M13 (Khalack et al., in preparation) suggests that 
the trends observed in HD 135485 for N and S, and possibly 
Fe, are real. These facts can be considered as further argu- 
ments in support of the efficiency of atomic diffusion in the 
atmosphere of HD 135485, and possibly in the atmospheres 
of BHB stars generally. 

The study of chemical stratification in the atmospheres 
of BHB stars is a recent field of research. Previously 
Bonifacio et al. (1995) reported the possibility of he- 
lium stratification in the atmosphere of Feige 86. The 
results obtained here for nitrogen and sulfur stratifica- 
tion in HD 135485 are the first report of stratification of 
metal species in a BHB star. An extensive search for ver- 
tical stratification for a large number of chemical species 
in the atmospheres of BHB stars with a range of T G ff 
would be very useful to place constraints on model atmo- 
spheres in which the atmospheric structure is calculated 
self-consistently with the stratification predicted by the dif- 
fusion phenomenon, such as those of Hui-Bon-Hoa, LcBlanc 
& Hauschildt (|2000|) . To obtain more precise results from 
observations pertaining to stratification of the elements, 
such models should be used since stratification can modify 
the atmospheric structure and thus the inferred stratifica- 
tion obtained by spectral analysis. For example, we found 
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that the nitrogen stratification for HD 135485 is quite sen- 
sitive to the T e g used for the underlying model atmosphere. 
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